Analysis and minimization of bending losses in discrete quantum networks 
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We study theoretically the transfer of quantum information along bends in two-dimensional dis- 
crete lattices. Our analysis shows that the fidelity of the transfer decreases considerably, as a result 
of interactions in the neighbourhood of the bend. It is also demonstrated that such losses can be 
controlled efficiently by the inclusion of a defect. The present results are of relevance to various 
physical implementations of quantum networks, where geometric imperfections with finite spatial 
extent may arise as a result of bending, residual stress, etc. 
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I. INTRODUCTION 

Studies on the faithful transfer of quantum information 
and the engineering of discrete quantum networks have 
been focused mainly on one-dimensional topologies, and 
the transfer of signals between the two ends of a quantum 
chain However, in analogy to conventional networks 
, the prospect of large-scale quantum information pro- 
cessing (QIP) and networking, irrespective of the physical 
platform, require efficient complex signal manipulations 
(such as routing, splitting, switching, etc), which are pos- 
sible only in higher-dimensional geometric arrangements. 
This necessity has motivated studies on state transfer 
in various 2D arrangements 0], most of which rely on 
nearest-neighbour (NN) Hamiltonians. 

Bends are expected to be at the core of any 2D con- 
figuration, yet their role on the fidelity of the transfer 
has not been investigated in the literature so far. We 
address this issue by investigating the transfer of quan- 
tum signals along a bent quantum chain whose operation 
is based on two different faithful-communication (FC) 
Hamiltonians with NN couplings [3, [f| , one of which has 
been employed in many of the 2D arrangements of Q. 
By construction, these Hamiltonians ensure the faithful 
communication through the unbent chain, whereas their 
performance in the presence of bends is not known, and 
will be discussed here. In many physical realizations, the 
coupling between two adjacent sites is directly related to 
their spatial separation (e.g., see [||-[8|). Hence, interac- 
tions beyond NNs are expected to get enhanced in the 
neighbourhood of bends, disturbing the communication. 
One way to circumvent such problems is to engineer new 
FC Hamiltonians for the bent chain taking into account 
interactions beyond NNs @ ; a rather tedious process for 
long chains. Here, we adopt another approach namely, 
the minimization of bending losses on the basis of the 
unperturbed NN Hamiltonians by means of minimal ex- 
ternal control (i.e., without elaborate sequences of pulses 
and measurements). To this end, a thorough analysis on 
the bending losses has been performed. 

In the following section we formulate the problem, 



whereas in Sec. Ill we analyse the bending losses and 
provide a way for their minimization. We conclude with 
a discussion in Sec. IV. 



II. FORMALISM 

The 2D arrangement under consideration pertains to 
N identical sites and is depicted in Fig. Q] Each site is 
associated with a qubit and the entire structure operates 
as a bent chain of qubits that interact according to a 
Hamiltonian of the form 



(la) 



where is the unperturbed Hamiltonian corresponding 
to the unbent wire i.e., to 6 = 0, and ^(6) is the per- 
turbation associated with the bending. We consider two 
protocols, which ensure the faithful transfer of informa- 
tion between the two ends of the unbent chain, and our 
task is to analyse their robustness against perturbations 
that stem from the bending of the chain. Both of the 
protocols pertain to centrosymmctric channels with NN 
interactions and the unperturbed FC Hamiltonian is of 
the form [<| 
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where a J is the creation operator for an excitation on the 

fijj is the 



jth site of the channel with energy Ej, and 
coupling between the sites with indices i and j. 

Protocol 1. The first protocol is characterized by Sj = 
e, fij,j+i = f^o y? 7^ 1, N — 1, and 0^2 = f2jv-i,iv = 
[5j. For a given N, the ratio fl/flo is chosen so that 
faithful transfer of information between the two ends of 
the chain occurs at time [loj ]. 

Protocol 2. The second protocol also pertains to reso- 
nant sites, while the couplings along the entire chain are 
engineered according to Hjj+i = Qq\/(N — j)j In 
contrast to protocol 1, this scheme promises ideally per- 
fect transfer at time = 7r/(2f2o). From now on, we 
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FIG. 1. (Color online) (a) A bent quantum chain of N sites, 
(b) A closeup of the bend. Green (solid) arrows denote NN 
interactions of strength Qi,j. Red (dashed) arrows denote 
interactions beyond NNs, with the strongest one (of strength 
g) corresponding to the first neighbours of the corner site a. 



also refer to a chain that operates in protocol i as chain 
(channel) i. 

In various implementations of quantum networks, the 
coupling constant VLij depends, among other parameters, 
on the inter-site separation fjj ; typically, it is expected to 
increase as we decrease r^j (e.g., see @t3)- The details of 
the dependence of f2jj on r^.j may vary from realization 
to realization, but also from site to site within the same 
network because of disorder and imperfections. As long 
as, however, such a type of disturbances are sufficiently 
weak (a necessary requirement for QIP) it is reasonable 
to assume that, to a good approximation the spatial de- 
pendence of riij is governed by a law that is universal 
for the entire network in a particular realization. 

Let a £ [2, N— 1] denote the index for the corner site at 
the bend. The separation between the non-neighbouring 
sites with indices a + k and a — k' (where k, k' = 1,2, . . .) 
decreases with increasing 6, while the separation between 
NNs is assumed to remain unaffected (see Fig. [T]). Hence, 
couplings between non-neighbouring sites are expected 
to be present for sufficiently sharp bends, and have to 
be taken into account. In an attempt to understand the 
effect of the bending on chains 1 and 2, throughout this 
work we focus on a regime of 6 where the perturbation is 
dominated by the coupling between the first-order neigh- 
bours of the corner site with indices a ± 1, while the 
effect of couplings between higher-order neighbours on 
the time scales of interest can be neglected. Hence, the 
perturbation in Eq. ([1]) can be chosen as 

f(6) = g(6){at_ 1 a a+1 + (lc) 

The regime of 9 for which such an assumption is valid, 
is intimately connected to the details of a given imple- 
mentation. In a more general context, the perturbation 
(ITcl) is expected to describe other types of undesirable ge- 
ometric spatially-localized imperfections associated e.g., 
with residual stresses [Tlj . in a particular realization of 
the chains. In cither case, the key point is that the dy- 
namics of the chains, are expected to be determined by 
the strength of the perturbation g, relative to the NN 
couplings around the corner, rather than the actual ori- 
gin of the perturbation. In the limit of weak perturba- 
tion (g -C ^a.Qipi) the performance of the bent chan- 
nel is expected to be rather close to the performance of 
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FIG. 2. (Color online) Performance of the bent chain f for 
various N. (a-d) The normalized probability of transfer Qi, 
for the bent chain without (bars), and with (x) corner defect. 
Corner sites: (a,b) a — 6; (c) a = f f ; (d) a = f3. (e) Time of 
transfer through the bent chain without (open symbols) and 
with (filed symbols) corner defect, (f) Optimal detunings of 
the defect corner site relative to the other sites 113 ■ 



the unbent one, whereas large deviations are expected 
for g > O a ct=F i. In view of the absence of related ex- 
periments, we keep the following analysis in a rather 
general theoretical framework, by introducing the ratio 
k = ff/O max [3, where O max = maxjOij} [H]. 

Following P, 0-0] , the chain is initially prepared in an 
eigenstate of M ', and the information to be transferred is 
encoded in the state of the first qubit. The Hamiltonian 
([1]) preserves the number of excitations and thus, the 
system is restricted to the one-excitation Hilbert space 
throughout its evolution. Various degrees of freedom 
that may be associated with the quantum state of the 
information carrier are assumed to be preserved on the 
time scales of interest, and thus the problem of the state 
transfer boils down to the transfer of the excitation. The 
computational basis can be chosen as where |j) 

is the state with one excitation on the jth site. In the 
Schrodinger picture the state of the system at any time 
t is given by = J2j A j li)> where \Aj(t)\ 2 is the 

probability for the excitation to occupy the jth site at 
time t. Initially, |^(0)) = |1), and the evolution of the 
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FIG. 4. (Color online) Performance of bent chains for different 
positions of the corner, and for fixed ft = 0.4. Different points 
for each TV are obtained by increasing a from 2 to [TV/2] with 
step 1 • (a,c) Protocol 1 and (b,d) Protocol 2. 
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terms of the ratios Q t = P l /Pj Q} and S t = Ti/T^ . 
of the results presented here pertain to a bend in the 
middle of the chains. Analogous observations and con- 
clusions hold for all a G [2,7V — 1], and thus we do not 
show related plots [HI]. 



FIG. 3. (Color online) As in Fig. [2] for chain 2. 



A. Analysis of bending losses 



amplitude Aj for protocol i and given a is governed by 
.dAj(t) 

i — — — - EjAj + S2j_ijAj_i + U jtj+ iAj +1 

+9^j,a+l-A a -l + gSj^a-iAa+i, (2) 

with j G [1,7V] and S m>n is the Kronecker's delta. The 
last two terms are associated with the bend, and they 
were not present in any of the previous investigations 
on the two protocols We are interested in the 

probability for the excitation to occupy the Nth site 
at time t when protocol i is used, which is given by 
Pi (t) = \Ajy(t)\ 2 . The corresponding probabilities for 

the unbent chains are denoted by i^ (0) (i). 

III. SIMULATIONS 

Equations ([2]) were solved numerically for various val- 
ues of a and n G [0, 1], keeping track of the first max- 
imum of Pi{t) and the corresponding time Tj, at which 
this is attained [l4[ ■ The evolution was restricted to times 

t G [0,T^ ], since the perturbation accelerates the trans- 
fer relative to the unbent channels (intuitively speak- 
ing, the coupling strength around the corner increases). 

Given that both of P 1 (0) and T 1 (0) are functions of N [|, 
for the sake of comparison our results are presented in 



According to the histograms of Figs. [2]and[3j for both 
protocols Qi w 1 for relatively weak perturbations (up to 
K w 0.2 ), and drops as we increase k. For a given proto- 
col, there do not seem to exist major differences between 
even and odd values of N, whereas the decrease becomes 
slightly faster as we decrease N. This is due to the spa- 
tial extent of the perturbation, which pertains to two 
sites around the corner site. The fraction of perturbed 
over unperturbed sites is thus getting smaller as we in- 
crease N, and the effect of the perturbation becomes less 
significant for a given value of k. A Gaussian fit turns 
out to be a rather good approximation for the estimated 
points Qi(n), and the ratio of the corresponding widths 
for the two protocols (associated with the decrease of Qi 
for increasing k), is estimated to about 1.08; a fact that 
shows how close the responses of the two protocols to the 
perturbation are. 

The effect of the bend on the time of transfer for var- 
ious iV is depicted in Figs, [^c) and[3Jc) (open symbols 
only). For both protocols the transfer is accelerated rel- 
ative to the unbent chains, with the acceleration being 
slightly more pronounced for chain 1. For a given pro- 
tocol, the acceleration becomes less pronounced as we 
increase TV; a behaviour that can be attributed again to 
the finite spatial extent of the perturbation. Moreover, 
the bend seems to affect mostly the fidelity and secondly 
the time of the transfer. 
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FIG. 5. (Color online) Difference between successive eigenen- 
ergies of the Hamiltonians for N = 25 and a — 12. (a) Pro- 
tocol 1, K = 0.5; (b) Protocol 2, k = 0.3. 



B. Minimization of bending losses 



to gain further insight into the role of the bends and 
the minimization of the associated losses, we have anal- 
ysed the spectrum of the Hamiltonians (one-excitation 
sector) for the various cases [l9|. In Fig. [3] we plot the 
separation between successive eigenvalues for the unper- 
turbed Hamiltonian as well as the perturbed Hamiltonian 
with and without optimisation. Clearly, the presence of 
the perturbation affects significantly the relative position 
of the eigenvalues, and this disturbance is responsible 
for the observed decrease of the probability of transfer. 
When the detuning of the corner site is optimized, how- 
ever, the initial distribution is restored to a large extent, 
minimizing thus the losses. Some deviations at the bor- 
ders are not of great importance since the overlap of the 
initial state [^(0)) with the corresponding eigenvectors is 
negligible (i.e., the overlap is peaked around the center). 



Various approaches for the minimization of the bend- 
ing losses have been employed. The most efficient we 
found pertains to the introduction of a defect by adjust- 
ing the energy of the corner site, while keeping all the 
other parameters the same. Let us denote by A Q , the 
detuning of the corner site relative to the other sites of 
the chains. This detuning is optimised so that the trans- 
fer from the first to the last site is maximized for a given 
strength of the perturbation k and for times t G [0, T^]. 
The optimal values of A a for the two protocols and var- 
ious N are depicted in Figs. |^f) and EI f) H3, while the 
corresponding values of Qi arc depicted with symbols ( x ) 
in Figs, [^a-d) and^a-d). Clearly, for given N, n and a, 
there is an optimal value of A a for which the probability 
of transfer is above 99% of the corresponding probability 
for the unbent chains fl8j . The optimal detuning turns 
out to be negative for all the tested parameters and in- 
creases (in absolute value) as we increase k. Surprisingly 
enough, for both protocols, the optimal values of A Q for 
various N do not differ substantially throughout the en- 
tire regime of k. Furthermore, to a good approximation, 
for both protocols and for all N, the behaviour of A a 
for k < 0.7 is linear, with estimated slopes —2.361 and 
—2.573 for protocols 1 and 2, respectively. In contrast to 
Qi, the effect of the optimization on Si is not so signifi- 
cant (filled symbols in Figs. [2{c) and(3Jc)); overall, the 
transfer is a bit faster relative to the perturbed chains 
without optimization (open symbols). 

According to Figs. B|a,b) the optimization is equally 
expensive for the two chains, in the sense that detun- 
ings of the same order are required in order to minimize 
the losses against the same disturbance K. It becomes, 
however, particularly expensive when the bend is close to 
the ends, and in this respect such arrangements should 
be avoided. The optimization seems to work efficiently 
irrespective of the position of the bend on the chains. 

A close inspection of Figs. HJe,f), 0He,f), and H{a,b) 
reveals that in general the optimal detunings are not so 
large to allow for adiabatic elimination of the corner site, 
and thus the reduction of the iV-site bent chain to an ef- 
fective (N— 1)— site chain with NN couplings only. Hence, 



IV. DISCUSSION 

We have discussed the transfer of quantum informa- 
tion along bent quantum chains that operate according 
to known FC Hamiltonians with NN interactions. Bends 
are at the core of various 2D configurations that have 
been discussed in the literature 0], yet their effects on 
the transfer have been neglected, and the related inves- 
tigations were focused on NN interaction Hamiltonians. 
Our analysis shows that the transfer is distorted signif- 
icantly by interactions beyond nearest-neighbours that 
stem from the bend. Nevertheless, the limited spatial 
extent of the perturbation, allowed us to minimize effi- 
ciently such losses by controlling the energy of the cor- 
ner site. Large-scale QIP requires reliable and efficient 
navigation of quantum signals in higher-dimensional net- 
works, where bends are expected to play a pivotal role. 
Our work sheds light on the role of such bends facilitating 
the engineering of reliable quantum networks in higher 
dimensions, including the 2D geometric arrangements of 

The present results have been obtained in a rather gen- 
eral theoretical framework, and are expected to be of rel- 
evance to various physical implementations of quantum 
networks pertaining e.g., to quantum dots, optical lat- 
tices or photonic lattices d, H, [2(| ■ In contrast to anal- 
ogous theoretical and experimental studies in the con- 
text of discrete soliton networks (see [2(| and references 
therein), our work pertains to linear networks whereas 
not all of the couplings between adjacent sites are the 
same. The versatility of photonic lattices, however, al- 
lows for engineering of various couplin g co nfigurations 
and geometric arrangements with bends [201 ] , so that our 
main observations can be confirmed experimentally in 
this context, with today's technology. 

In the weak-coupling regime the coupling between the 
waveguides with indices j and (J + 1) depends exponen- 
tially on their separation dj j+i i.e., we have Cjj+i = 
rj exp (— £dj.j + i) , where rj, £ are open parameters to be 
determined by fitting to related experimental data for 
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a particular setup |8j. Given this exponential law, one 
can readily show that a particular coupling distribution 
is obtained in practise if the distance between 
successive waveguides is chosen according to djj+i = 
— In m^fijj+i) As typical values for the open pa- 
rameters at the wavelength A = 633 nm, we may consider 
r\ ~ 19.5 cm^ 1 and £ = 0.152 pm" 1 . Hence, the separa- 
tions required e.g., for the realization of the coupling dis- 
tribution of protocol 2 in a lattice of N = 9 waveguides 
of length L = 10 cm, range from about 21.86 to 24.88 
fim. Such separations are well within reach of current 
technology used in fabricating photonic lattices @. In 
view of the above exponential law, one can also readily 
show that NN couplings are almost two orders of mag- 
nitude larger than the couplings beyond nearest neigh- 
bours. One can thus fabricate higher-dimensional array 
configurations with NN couplings, but in the neighbour- 
hood of bends interactions beyond NN set in locally ■ 
In the context of waveguides, the detuning of the cor- 
ner site relative to the others can be achieved by altering 
the core refractive index or the core radius of the corner 
waveguide along the lines of [2(| ■ Such an adjustment has 



been demonstrated and studied experimentally for bent 
photonic arrays recently pll ] . As mentioned before, how- 
ever, these studies pertain to discrete soliton networks, 
whereas our work paves the way for manufacturing reli- 
able higher-dimensional linear networks with engineered 
couplings, by minimizing losses associated with bends. 

In closing, we would like to emphasize that in princi- 
ple there are infinitely may state-transfer Hamiltonians 
p}, and the present analysis can be performed for any 
one of these Hamiltonians as well. Moreover for a given 
Hamiltonian, one may explore many different schemes for 
minimizing bending losses. Hence, the comparison of dif- 
ferent Hamiltonians and different minimization schemes 
can be only a long-term project that goes clearly beyond 
the scope of the present manuscript. 
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